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In high energy nuclear collisions, the Glauber model is commonly used to evaluate J/ip sup- 
pression due to inelastic interaction with colorless bound nucleons. This requires an effective 
^ ' value for the J/?/>- nucleon absorption cross section which is larger than theoretically expected. 

On the other hand, multiple nucleon-nucleon collisions mediated by color exchange interac- 
tions, excite their color degrees of freedom. We investigate the importance of this effect and 
conclude that the related corrections are important to explain the effective value extrapolated 
from experiment. 
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In a nucleus-nucleus (AB) collision JM suppression is mostly treated within the Glauber 
model, with a constant effective cross section for absorption on bound, non-interacting nucleons. 
Here we intend to study and to improve the understanding of the physical phenomena behind 
early stage absorption. Following a known interpretation of the dynamics of nucleon-nucleon 
(NN) collisions, we ascribe the main contribution to the inelastic cross section to color exchange 
processes. Via color exchange, the 3-quark systems constituting the colliding nucleons become 
color octet states. Color exchange may also be accompanied by eluon radiation. While the 
contribution of radiated gluons to J/ip absorption was treated in El, here we focus on the first 
part. Having realized that color-singlet nucleons are found in color octet states after an elemen- 
tary collision, one can expect that repeated scatterings lead to multiple color exchanges. One 
therefore has the possibility to excite the color degrees of freedom of the 3-quark system also 
^ . to the decuplet state. Higher color representations can be obtained when radiated gluons are 

also present, but here we restrict our study on the color dynamics of the constituent quarks of 
a nucleon. 

To calculate the value of the inelastic cross section for f ona nucleon which has undergone 
multiple scattering, we use the two-gluon exchange model of Low and NussinovH Although 
oversimplified, the model incorporates important properties of high-energy hadronic cross sec- 
tions such as (approximate) energy independence and scaling of the cross sections according to 
the mean square radius of hadrons. Indicating with the effective charmonium state, composed 
by a mixture of J/ip, ip' an d X> one can write at once the absorption cross section, obtaining 
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The first factor in the integrand is the squared gluon propagator, the second contains the phase 
shift due to the momentum transfer k and is responsible, by the way, for color transparency, 
while the last factor is the interesting part of the present calculation. It is the product of 
two vertices between three quarks and a gluon, averaged with the wave function of the colored 
nucleon. The vertex is T a (k) = 16/3 a s [ A"^) exp( i k y% ) + A" 2 ) exp( ikfa) + X%\ exp( i ky^ )]. 
In order to average the product of the two vertices, one needs to specify the structure of the 



nucleon wave function. The color part of the wave function is the relevant one, while the spatial 
part is chosen to be equal for all color multiplets. The factor given by the angular brackets in 
eq. (|l]) can be evaluated and leads to 

(N n \r a r^ a \N n ) = (N\J2 l-CneMikVij^N). (2) 

i<j 

The effect of having different color states is all contained in the coefficients c n , listed in the 
second column of Tab. [I]. The values of c n found, for n = 8, 10, lead to larger absorption cross 
sections as compared to the singlet case. To see this it is useful to define the two-quark form 
factor F^(k 2 ) = (N\ex.p(ikxij)\N) for the nucleon. Since, within our approximation, the 
latter does not depend on the colour multiplet, we drop the label n. Then we can rewrite eq. ([l]) 
as a$fc = (\&| a n (r) where 

*n(r) = f a 2 s j d 2 k ^ 2)2 [1 - exp(ifcF)] [1 - c n F 2 N q (k 2 )] , (3) 

exhibiting a similar structure to that of the usual dipole-nucleon cross section, but with the 
important difference carried by the coefficients c n . Notice that with c n ^ 1 and \x = the 
integral in eq. (||) is infrared divergent. Therefore, the cut-off /j,, of the order of the inverse 
confinement radius, is important. With /j, ^= it holds that <r n (r) — > for r — > 0, so that one 
still preserves color transparency, even with a colored nucleon. This is indeed a property of the 
singlet meson. Here we are not primarily interested in the absolute values of but rather 

on its dependence on the color state n with respect to the singlet value. Yet we estimated the 
absolute values choosing a s = 0.6 and [i = 140 MeV. Within a harmonic oscillator model of the 
nucleon, one finds that F^(k 2 ) = F^J n (3k 2 ), the electromagnetic form factor. We assume that 
this property is general and take the dipole form F^iQ 2 ) ~ (A 2 + Q 2 )~ 2 with the parameter 
A 2 = 12/ R 2 and R p = 0.8 fm. With the above definitions, the generalized dipole-nucleon cross 
section can be computed analytically and the results of the calculation are illustrated in Fig. ??, 
showing the modified dipole-nucleon cross section for the three color states of the nucleon. It 
largely increases when the dipole scatters off a colored object. One can get more insight in the 
effect of the different color states rewriting eq. (|3|) in the form 

a n (r)=a(r)[l + (l-c n )A(r)\ , (4) 

where a(r) is the usual dipole cross section on a color singlet nucleon and the function A(r) is 
found to vary appreciably with r. Its functional dependence is given in the small box in Fig. 
??. To calculate the inelastic cross section, we must average a n with the meson wave function. 
Therefore, the dependence of A on r shows that the effect of having colored nucleons is stronger 
for larger size mesons. Recent calculations i, based on realistic phenomenology for the dipole 
cross section and charmonium wave function, gave values for cross sections of various charmonia 
on a nucleon. At center of mass energy = 10 GeV, it was found that (Jjj^ n = 3.6 mb, 
&ip'N = 12.2 mb and o x n = 9.1 mb. With a composition of ^ of 52-60% due to J/tp, 8-10% due 
to V' and 32-40% due to \i a weighted value of 5.8 mb was obtained. In the present case we 
evaluate with a Gaussian wave function having root mean square transverse separation 



Table 1: Results of the calculation of the meson-colored nucleon cross sections. 
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-1/4 


7.8 


1.35 


Decuplet 


{10} 


1/2 


9.9 


1.7 



(v 2 ) 1 / 2 = 0.25 fm, in order to account for the different states implicitly included in and 
therefore adjust the singlet cross section to the quoted 5.8 mb. The values found are given 
in the third column of Tab. [l|. In order to minimize the model dependence, we factored out 
from the colored cross sections the amount corresponding to the singlet case, therefore writing 
a $N n = a< ^N The values for 8 n are given in the last column of Tab. [j]. They constitute the 
important result of this calculation , showing that ^ scatters with an color octet nucleon with a 
35 % larger cross section with respect to a singlet nucleon, while with the decuplet the increase 
is even of 70%. The large values found now require a careful analysis to establish the validity 
of this newly proposed absorption mechanism. 

Since the absorption pattern observed in pA collisions is usually considered to be the baseline 
for the AB case, it is necessary to examine it in some detail. Before we analyze the data for 
charmonium production in pA collisions, we emphasize the importance of charmonium formation 
time ? . Within a simplified version of a two-channel model □, the N absorption cross section 
evolves in the eigentime r of ^ as a^(r) = o % ^ N + (a^ — a % $ N ) cos(t/t/), where t/ = (m^i — 
m Jijj)~ l = 0.3 fm is the charmonium formation time, (7^ N is the inelastic cross section and 
is the cross section of the so callecL "pre- meson" on a nucleon. As mentioned in the last 
section, calculations give cr^ N = 5.8 mb@. Unfortunately has not yet been calculated. We 
therefore use the value 2.7 mb fromi. 

One can then write the expression for ^ production in pA collisions, as 
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We want to compare this expression with the conventional one obtained in the Glauber model, 
containing only an effective absorption cross section, therefore making the replacement <t|^(t) — > 
cj^y. For A = 208 and p^ = 50 GeV, using a realistic profile for Pb we obtain cr^y = 3.8 mb. 
This number has to be interpreted as the value of the effective absorption cross section of ^ on 
a color singlet nucleon, including formation time effects. To obtain the corresponding value for 
the case when \E' scatters on a color octet nucleon we have to enhance the singlet result by 35%, 
obtaining 5.1 mb. 

We have then re-analyzed the existing data taken by 3fixed target experiments: NA3 at 
CERNQ, NA38/51 i at CERN and E866 at FERMILAB I. In order remove the effects on 
absorption due to different \£ momenta, we impose the same kinematic conditions, considering 
^'s of 50 GeV momentum in the laboratory frame for the three experiments. We then assume 
that the absorption mechanism is parameterized by the effective cross section (T^t. We made 
three separate fits including quoted statistical and systematic errors and extracted the effective 
absorption cross section by minimizing x 2 - The result of these fits is presented in the second 
column of Tab. ^. The mean value of 3.5 db 0.6 mb seems to favor the singlet mechanism, for 
which a cross section of 3.8 mb is theoretically expected. 



Table 2: Results of the fits to pA data for the various experiments considered. See text for details. 
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offi [mb] 


Glauber fit 


Gluon corr. 


NA3 


4.5 ± 1.7 


7.3 ± 1.9 


NA38/51 


7.1 ± 1.6 


10 ± 1.8 


E866 


2.7 ± 0.7 


4.7 ± 0.7 


Averages 


3.5 ± 0.6 


5.7 ± 0.6 


Theory 


3.8 {1} 5.1 {8} 



On the other hand, the assumption that ^ is produced on a bound nucleon with the same 
cross section as on a free one, can be an oversimplification of the physics under study. In a nuclear 
environment parton distributions are substantially modified. In particular, the gluon distribution 
in nuclei exhibits anti-shadowing effects at x B j — 0.1. We expect that charmonium production on 
a nucleus is enhanced by a factor R g A (x2, Q 2 ) = gA(x2, Q 2 )/d(x2, Q 2 )- We estimate these effects 
by rescaling the ^ cross section with the A-dependent correction R 9 A = 0.032 log(vl) + 1.006, 
extracted fromB. We fitted the data once more, obtaining the results listed in the third column 
of Tab. One sees that the new average value of 5.7 ± 0.6 mb now seems to favor the octet 
mechanism, for which we estimated a value of 5.1 mb. One sees that the situation is, in all, 
quite uncertain. 

The above analysis assumes that there is no dependence of nuclear effects on the proton 
energy if the charmonium energy is the same in all experiments. However, there might still be 
some dependence on the proton energy. According toS, energy loss of projectile partons in the 
nuclear medium substantially reduces the production rate of charmonia at large xf and is still 
rather important at small xf- In order to compensate for energy loss, the gluon distribution 
must be evaluated at a higher value of x±, therefore suppressing \& production. This will in turn 
reduce the final state absorption on nucleons in a fit to the data. To account for this effect, 
further study is required. 

To summarize, the data seem to suggest that ^ is absorbed by colored nucleons in pA colli- 
sions. The picture is though not yet conclusive due to the large experimental errors. Neverthe- 
less, we follow our theoretical prediction and examine the potential role of the color mechanism 
in the more complex case of AB collisions. Preparing the ground for the forthcoming discussion, 
it is useful to recall the conventional expression within the Glauber model for the $ production 
cross section at impact parameter b. One has 



— 4^(6) = d 2 s dzAdzB cr^ N Pa(za,s) PB(zB,b - s) 
d 2 b J 

1 - a% T A (z A , s)/A] A ^ [l - af s N T B (z B ,b- s)/b] *~\ (6) 
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Here we neglect the effects of any additional suppression factors which go beyond conventional 
nuclear effects, that is final state interactions with the produced matter (QGP, hadron gas, etc.), 
since we are not primarily interested in these contributions. On the other hand, we focus on the 
reinterpretation of the terms describing nuclear effects discussing color exchange in the context 
of multiple scattering, then applying the obtained results to the calculation of production in 
AB collisions. 

Having completed the computation of the \P cross sections with colored nucleons, we need 
to look at how these color objects develop in a collision between nuclei. Consider a nucleon of 
one of the nuclei colliding with a row of nucleons of the other nucleus. It undergoes multiple 
scattering, with repeated color exchanges. It is possible to calculate the probability P n (z,b) for 
a nucleon in color state n = 1,8,10 after having traveled through the row up to z at impact 



parameter b. This can be done in a very e 
color density matrix of the 3-quark system 



egant way, by solving an evolution equation for the 
L9. One obtains 



, 1 + 20 F(z,b)+ 6 G(z,b) {1} 
Pn(z,b) = —x\ 16-40 F(z,b) + 24 G(z,b) {8} , (7) 



10 + 20 F(z,b) -30 G(z,b) {10} 



where F(z,b) = exp(— g cr l ff N T(z,b)) and G(z,b) = exp( — 7 cr l $ N T(z,b)). One notices several 
properties of the probabilities. First of all, the limit T — > implies F,G —* 1. This means that 
P\ —* \ and Ps,P\o ~~ > 0. In other words the nucleon, before entering the nucleus, is in singlet 



state as it should be. Moreover, P± + Ps + Pio = 1 for any T, therefore probability is always 
conserved. Finally, if the nucleus is large enough one has T » 1 and F, G <C 1. This implies 
that if sufficient scatterings have taken place, the color probabilities reach the statistical limit, 
given by the first coefficients of eqs. (|7|). It turns out that soon after the nucleon has penetrated 
the nuclear profile, the statistical limit is essentially reached in which about 16/27 ~ 2/3 of 
the nucleons are in octet state while 10/27 ~ 1/3 are in decuplet. The amount of singlet is 
negligible. 

With the calculated probabilities P n and the previously evaluated color cross sections er^y , 
we can now construct the cross section for ^ when it scatters with a nucleon, taking into account 
all its color states, and modify eq. (^) by making the replacement 

dz A E B (z B , b B ) Pa(z' a , b A ) , where z B = z B - (z' A - z A ) - • (8) 

-oo I — Via 

The position-dependent effective cross section 

X B (z B ,b B ) =E ot** Pn(*B,b B ), (9) 
n 

represents absorption by nucleus B, taking into account the multiple scattering of its nucleons 
with nucleus A. Analogous expressions hold for A. 

The modified version of eq. @ can now be used to evaluate total cross sections, with a 
further integration in impact parameter. We compute^alues to compare with available data 
measured by the NA38/50/51 collaboration at the SPS". We use realistic nuclear profiles and 
scale the center of mass energy according to the parameterization a^ N ~ (1 — rrij/^/y/s) 12 . 
The results are compared with the data in Fig. |]. The dashed curve corresponds to a standard 
Glauber model calculation with the singlet absorption cross section (T^ft = 3.8 mb, as evaluated 
in Sect. 3, while the thick curve refers to the improved calculation which takes into account color 
excitations, studied in Sect. 4. In order to emphasize the difference between the two absorption 
mechanisms, we neglect here anti-shadowing effects. For the pA calculation we increased the 
value of cr'yxi by 35%, therefore used (T^jy = 5.1 mb. This is the extreme situation, when \I> 
always encounters color octet nucleons. For the AB case we used eq. 

Different from the Glauber model result, our improved one provides a better agreement with 
the data, except for the Pb+Pb point. In other words, the proposed absorption mechanism due 
to color dynamics is found to be an important contribution to nuclear effects in ^ production. 
It must be accounted for, together with other mechanisms of suppression which we here put 
aside. 

To summarize, we studied the dynamics underlying ^ suppression due to nuclear effects, 
employing the idea that color exchange processes in multiple NN collisions may play a significant 
role. We calculated inelastic cross sections for ^ scattering on colored nucleons and found an 
increase of about 35 % when the nucleon is in a color octet state as compared with the singlet 
case, while the enhancement factor for the color decuplet is about 70 %. 

Making use of available calculations of *$>N cross sections, we included formation time effects 
and obtained the effective value for absorption on singlet nucleons = 3.8 mb. Scaling 

this value with a 35 % increase we obtained cf, 6 ^ = 5.1 mb and with a 70 % increase c^io = ^ 
mb. We then performed a fit to available data on ^ production cross sections in pA collisions. 
The result of the fit, including gluon anti-shadowing corrections, is that the effective absorption 
cross section for ^ is (J^t = 5.7±0.6 mb, therefore suggesting that in pA collisions, the produced 
charmonium interacts with color octet nucleons. The conclusion is, though, rather weak because 
the experimental points are quite scattered. 

Total production cross sections for were then calculated and compared in Fig. |l| with 
available data for AB collisions. Compared to the Glauber model result, the agreement is 
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Figure 1: Left: Dipole cross section for different color states as a function of the quark separation r. Box: The 
quantity A(r) which governs the dependence on color. Right: Calculations of total J/ip production cross sections 
in AB collisions and comparison with the NA38/50/51 data. The curves, drawn to guide the eye, correspond to 
Glauber absorption (dashed) and to absorption by colored nucleons (full). 



improved and we conclude that effects due to color exchange processes are substantial and 
provide the value cr^fe = 5.5 mb for the effective absorption cross section on nucleons in AB 
collisions (See eq. (S)). 

In conclusion, we have completed a study of nuclear effects in charmonium production in 
AB collisions, providing a better baseline for future calculations. These amount to an improved 
treatment of gluon bremsstrahlung, including non-linearities arising because of gluon cascades, 
to the calculation of effects caused by gluon energy loss and to the evaluation of -E^-dependent 
cross sections. This will hopefully provide the starting point from which to address more precisely 
the nature of anomalous suppression in Pb+Pb collisions. 
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